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We discuss a fundamental effect of the interaction-induced decoherence of the electron wave 
function in disordered metals. In the first part of the paper we consider a simple model of a 
quantum particle interacting with a bath of harmonic oscillators and analyze the physical origin of 
the effect. This exactly solvable model also allows to understand why the arguments against the 
existence of the effect at low temperatures fail. The second part of the paper is devoted to a rigorous 
analysis of quantum decoherence in disordered metals. We also discuss the relation of our results 
to the recent experiments on GaAs structures. The existence of a finite quantum decoherence rate 
at low T implies that low dimensional disordered metals with generic parameters do not become 
insulators even at T = 0. 

I. INTRODUCTION 

The concept of quantum coherence is one of the most fundamental in quantum mechanics. According to the general 
principles quantum coherence of the wave function cannot be destroyed due to elastic interaction with a static external 
potential. On the other hand, inelastic scattering processes may (and in general do) destroy the phase coherence of the 
wave function. For electrons in a metal such processes (e.g. inelastic electron-electron and electron-phonon scattering) 
are important at sufficiently high temperatures. But as the temperature T is lowered inelastic processes become less 
intensive and the corresponding inelastic scattering time Tj for an electron in equilibrium tends to infinity at T — > 0. 
Hence, one could consider natural that interaction may cause dephasing in a metal at nonzero T, but not at T = 0. 

Surprizingly, it was found in many experiments with disordered metals (see jl] and references therein) that the 
effective dephasing time t v for the electron wave function saturates at the level t v ~ 10 -1 -j- 10 ns and does not 
depend on temperature below T ~ IK. These findings are in a clear contradiction with the above point of view 
according to which the time t v should increase to infinity as T approaches zero. Scattering on magnetic impurities, 
heating and the external noise were suggested as possible reasons for saturation of t v observed in earlier experiments. 
All these reasons were convincingly ruled out by the authors JjJ] , and it was argued that it is the zero point motion of 
electrons in a disordered conductor that causes dephasing at T = 0. 

A theory of the effect of quantum decoherence at low temperatures was suggested by the present authors g. It 
was demonstrated that at sufficiently low T the high frequency quantum noise of the effective electronic environment 
is responsible for the effect of quantum decoherence. It was predicted that at low T the decoherence rate l/r„ for Id 
metallic systems saturates at the level 

1/T V ~ e 2 t>F/7T(7l, (1) 

where o~\ is the wire conductance per unit length, vp is the Fermi velocity. The result ([!]) is in a very good agreement 
with the experimental findings as well as with previous experimental results. 

The predictions B inevitably lead to another fundamental conclusion: low dimensional metals do not become 
insulators even atT — 0. Indeed, e.g. in the case of Id systems the localization length is Zi oc ~ N c hl, while the result 
(|l|) yields the effective decoherence length L v ~ l\/N c h, where I is the elastic mean free path and iV c h is the effective 
number of conducting channels. It is obvious that for N c ^ 3> 1 (which is always the case in metallic wires) L<p is 
parametrically smaller than ?i oc . This implies that due to interaction with other electrons the phase coherence of the 
electron wave function is destroyed faster than the electron can get localized, i.e. localization never takes place under 
the above conditions. 

Since the results |l],|2| imply a necessity to strongly reconsider the commonly adopted point of view on the role of 
interaction in disordered metals at low T, it is not surprizing that the subject still remains controversial. Moreover, it 
is sometimes argued that quantum decoherence at T = would contradict to general principles of quantum mechanics. 
What are the arguments against the quantum noise in a disordered metal as a reason for quantum decoherence of 
electrons? Usually it is argued on a general level, adopting a much simpler quantum mechanical model in order to 
illustrate the main statement. On one hand this makes sense indeed: the nature of a fundamental quantum mechanical 
effect (or the absence of it) can and should be understood without unnecessary complications. On the other hand, one 
has to make sure that all significant features which yield the effect in a more complicated problem are still present in 
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a simplified one. There is an obvious danger in making too many simplifications: one can easily throw out the baby 
with the bathwater. 

Consider a quantum particle colliding with a harmonic oscillator with the frequency lu. Assume that the interaction 
potential is short range and before the collision the particle was in the state with the energy smaller than huj, while 
the oscillator was in the ground state. The general features of this scattering problem are well known: after interaction 
the oscillator remains in the ground state (the particle does not have enough energy to excite it) , the particle energy 
also remains the same as before interaction, the particle wave function - although in general changes as a result of 
interaction - stays fully coherent. Thus no energy exchange between the particle and the oscillator has occured, and 
no phase coherence has been lost due to interaction. What remains is to understand if the above consideration is 
relevant to the problem in question. 

At the first sight it appears to be relevant indeed. Although in a metal the electron interacts with other electrons 
instead of oscillators (let us ignore the electron-phonon interaction for simplicity), at T = all states below/above 
the Fermi energy should be occupied/empty, and no energy can be transferred: the electron cannot be scattered into 
any lower energy state due to the Pauli principle as well as to any higher energy state because no energy can be 
extracted from other electrons due to the same reason. It follows immediately, that at T — > the phase coherence of 
the electron wave function in a metal cannot be lost due to the very same reason as in the above example with the 
particle and the oscillator: there exists no inelastic scattering. 

This consideration turns out to be too naive, however. It may be sufficient if the infinite time limit can be taken 
already at the very first step of the calculation as e.g. within the standard formulation of the scattering problem. In 
this case the result will be proportional to the delta-function of the energy difference between the initial and the final 
states ensuring the energy conservation. However, at any finite time the energy-time uncertainty principle should be 
taken into account. The energy uncertainty can be not necessarily small, especially if the particle interacts with an 
infinite number of oscillators. This is crucially important for the effect in question. On top of that our problem is 
more complicated because we are dealing with an interacting quantum many body system rather than the scattering 
problem for two particles with well defined in- and out-states. The electron energy fluctuates as a result of interaction 
with other electrons, only the total energy of all electrons plus the interaction energy is conserved. One may introduce 
the effective oscillators also in this case, they are just the modes of the electromagnetic field. The interaction between 
the electron and one such mode can be naively described as H mt — ei?cos(kx), where E is the amplitude of the 
electric field, k is the wave vector and x is the electron coordinate. The interaction potential remains finite for any 
x, i.e. the electron always interacts with the oscillator and the scattering problem cannot be formulated. Actually 
the problem is even more complicated because the electron interacts with an infinite number of such oscillators at 
the same time. It is also important that in the presence of interaction the energy exchange between the low energy 
particle and the high energy oscillator is possible without excitation of the latter simply because its lowest energy level 
acquires a finite width. This effect depends on the strength of interaction, but it always exists because the interaction 
is never "turned off" . 

One might argue that this simply means that electrons are "bad" particles in the presence of interaction and one 
should rather define "better behaving" quasiparticles and calculate all measurable quantities in their terms. This does 
not always help because of at least two reasons. One is that quasiparticles (if they exist) are usually defined within 
an approximate procedure. Although the approximation may work well for calculation of certain physical quantities, 
it may fail for other quantities. Another reason is that the transformation of the basis (even if it can be done exactly) 
may not be convenient if the measurements are done only with a "bad" particle. In this case calculation with "good" 
quasiparticles can be by far more complicated, and it is much simpler to "get rid" of all but one "interesting" degree 
of freedom at an early stage of calculation by tracing them out in the full density matrix. This is the key idea of the 
Feynman- Vernon theory of the influence functionals (^^] developed further by Caldeira and Leggett j|] and Schmid 
Q in application to an infinite bath of harmonic oscillators with the ohmic spectrum. 

Superconducting Josephson junctions may serve as an example of a fermionic system where the same ideas have 
been worked out @||. In this case it is possible to exactly integrate out all electron degrees of freedom and describe 
the system dynamics in terms of only one collective variable - the Josephson phase ip. Of course, (p is known to 
behave "badly" (its quantum dynamics is incoherent in almost all cases, see e.g. ^-0) but namely this variable is 
of interest in Josephson junctions and SQUIDs just because the junction current and voltage operators as well as the 
flux operator in SQUIDs are defined in terms of (p. In this case "better behaving" quasiparticles are simply irrelevant. 
Therefore it appears to be unreasonable to argue in favour of "coherent" quasiparticles if all measurements are done 
only with "incoherent" variables. We believe that a similar situation is encountered if one calculates the conductance 
of a disordered metal. Actually in this case it is even more complicated because it is not clear if "good" quasiparticles 
can be introduced at all. 

The paper is organized as follows. In Section 2 we will demonstrate how the fundamental effect of quantum 
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decoherence can be derived within a simple model of a quantum mechanical particle interacting with an environment 
consisting of a collection of harmonic oscillators. It is remarkable that already this simple model captures all significant 
features of the effect and allows to understand why the arguments against its existence fail. We also establish the 
relation between our results and the standard perturbative treatment of a scattering problem. Section 3 is devoted 
to a rigorous microscopic analysis of quantum decoherence in a disordered metal. Discussion of recent experiments 
and comparison with our theory are presented in Section 4. The main conclusions are outlined in Section 5. 

II. FEYNMAN- VERNON THEORY AND QUANTUM DECOHERENCE 

A. Influence Functionals 

Consider a quantum-mechanical particle (or a more general system) characterized by a coodinate q(t) and the 
action So[q(t)]. Assume that this particle interacts with another quantum system described by a coordinate Q(t) 
and the action S env (Q). We will call the latter quantum system "environment". The total action for the system 
"particle+environment" has the form 

S[q,Q] = S [q] + S cnv [Q(t)} +S int [q(t),Q(t)}, (2) 

where S- m t describes interaction between q and Q. 

Suppose we are interested in the probability W for the particle q to have the coordinate q = qf at a time t provided 
at t = it had the coordinate q — qi. A general and elegant way to solve this problem can be formulated within the 
Feynman- Vernon theory of the so-called influence functionals. An extensive discussion of this theory can be found in 
Refs. [HJJ]. Here we only repeat the key steps. 

By definition the probability W is given by the square of a transition amplitude W — \K(qf, t; qi, 0)| 2 . In the absence 
of interaction (C = 0) the oscillator coordinate Q does not enter into the expression for K. However for C ^ there 
appears an additional force acting on a particle q. This force depends on Q which is itself a quantum variable. Thus 
we cannot anymore restrict ourselves to the dynamics of q, but rather should deal with the total Hamiltonian for a 
system il q + Q" or, equivalently, with a complete action (|J). Of course, the probability W is again equal to the square 
of the corresponding transition amplitude K which now depends on both q and Q. It is important, however, that we 
are not interested in the final state of the oscillator and do not make any measurements of Q. Therefore we should 
sum over all possible final states of Q and the formula for the probability W takes the form 

W = J2 K (lf> Qf> *i Qi,0)K*(q f ,Q f , t; e», Q 4 , 0). (3) 
Qf 

This is an important point. Each of the terms in the sum (||) represents a probability for a system to come into the 
state (qf,Qf). Since the subsystem Q can be in any final state Qf, in order to find the total probability W we have 
to add all these probabilities together. 

One can slightly generalize the problem and describe the evolution of the density matrix of the system p(q, q') from 
some initial to some final state. This evolution is described by the equation 

p(Qf,<l'f) = J dQidqiJ{qf,q'f,t;qi,q' i ,0)p i (q i ,q' i ), (4) 

where Pi(qi, q^) is the initial density matrix of the particle q. For the sake of simplicity in what follows we will assume 
that there is no interaction between q and Q before t — 0. Then the total initial density matrix can be factorized as 
Pi( a i,<li)Pi{Qi7 Q'i)- Here the kernel J is again given by the product of two amplitudes. It is convenient to represent 
this product in terms of a double path integral 

rif ri's 

J= D qi Dq 2 exp(iSo[qi(t)}-iS [q 2 (t)])F[ qi (t),q 2 (t)], (5) 

J qi J q'. 

where F is the influence functional which describes the total effect of the subsystem Q on the particle q. The functional 
F in turn can be represented in terms of a double path integral 
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F[ qi (t),q 2 (t)]=J2 J dQ, J dQ- DQi J Qf BQ 2 p(Q l ,Q' i ) 



rQ f fQ'f 

/ 

exp(iS'onv[(3l] -i5"cnv[Q2] + iSint [<?1 , Ql] - iSint[<?2,Q2])- (6) 



Here again (we cite from Q|) "X)/ j us ^ means that at some final time tf after we are no longer interested in the 
interaction we must take Qf = Q'f and integrate over all Q/", i.e. 

JdQ f J dQ' f 8(Q f - Q' f ). (7) 

This completes the general analysis of Feynman and Vernon. Its important advantage is that no approximations have 
been done so far, i.e. the above formulas are exact. One more citation from j| is in order: "F contains the entire 
effect of the environment including the change in behavior of the environment resulting from reaction with q. In the 
classical analogue, F would correspond to knowing not only what the force is as a function of time, but also what it 
would be for every possible motion q(t) of the object. The force for a given environmental system depends in general 
on the motion of q, of course, since the environmental system is affected by interaction with the system of interest 
q n . Thus all changes of the state of the environment resulting from the interaction with a dynamical variable q are 
automatically taken into account within the above formalism. 

In order to proceed we first assume that the environment consists of a single harmonic oscillator Q which has a 
unit mass and a frequency u>. For simplicity, the interaction is chosen bilinear with respect to both the particle and 
the oscillator coordinates q and Q, so that the total action for the system "particle+oscillator" has the form 

S[q, Q] = S [q] +Jdt Qq 2 - y<2 2 + CqQ^j , (8) 

where C is a constant which governs the strength of interaction. We will also assume that the oscillator is initially 
kept at a temperature T, i.e. the probability to occupy the state k is Wk = cxp(— kw/T)(l — exp(— uj/T)). In the limit 
T<Cw the initial state of the oscillator is its ground state k = 0. 

For this model it is a matter of a simple integration over the Q-variables to obtain the exact expression for the 
influence functional F (||). One finds HQ: 

F[qt,Q2] = exp(-iS R [q 1 ,q 2 ] - Si [qi,q 2 ])- (9) 

Defining q + = (q\ + q 2 )/2 and g_ = q\ — q 2 we have 

C 2 

Sr = — dtid« 2 g_(ti)sin(w(ti-i 2 ))g + (* 2 ), (10) 



C 2 / u \ r* 



Si = — coth {—) j dt 1 dt 2 q-(t 1 )cos(( J j(t 1 -t 2 ))q-(t2). (11) 



Let us point out that Si > for all trajectories q(t). Eqs. (|6[|ll|) summarize the complete effect of interaction with 
the oscillator Q on quantum dynamics of a particle q. 



B. Free particle interacting with oscillators 

Let us come back to the probability W (J3j) or, more generally, to the kernel J (||). We will consider two simple 
examples. The first example is a freely propagating quantum particle with a mass m. In the absence of interaction 
(C = 0) we have F = 1, and the double path integral (|5|) decouples into the product of two single integrals. For 
Qi/f = l'i/f cacn °f t nem is dominated by the same classical path q(t') = qt + (qf — qC)t' jt and (cf. e.g. Q) 

K(q f ,t;q t ,0) = y/m/2nitex.p(iS [q\), S [q] = m(q f - q t ) 2 )/t. (12) 

The probability W — KK* does not depend on the phase \S$ of each of the amplitudes, these phases enter with the 
opposite signs and cancel. 
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Now let us turn on the interaction (C / 0). If - just for the sake of simplicity - one treats this interaction 
perturbatively, one immediately observes that the double path integral (||) is again dominated by the same classical 
path qi(t') = q2{t') = q(t')- For any q\{t') = q2{t') we have Sr = Si = and therefore F[q,q] = 1 like in the 
absence of interaction. One can also establish the general form of the kernel J (0) . It is relatively complex and is not 
presented here. More interesting situation emerges if we modify our model and consider our particle q interacting with 
N oscillators with frequencies uj n . The corresponding generalization is trivial: one should just substitute w — ► uj n 



in (22) and carry out the summation over all 1 < n < N. If one sends N to infinity and assumes a continuous 
distribution of the oscillator frequencies H|(| : 

7rC 2 

V"- — [5{u n - u) - 8(u n + u)]=r}u, \w\<w c , (13) 



(uj c defines the high frequency cutoff) one arrives at the influence functional of the form l|) where one should 
substitute 

C \ n T c "duo 
w J- Uc 2tt 

The problem defined by eqs. (P- pr|Jl^ ) is gaussian and can be solved exactly, see e.g. jyj. Performing a straight- 
forward gaussian integration over q one arrives at the the exact expression for the kernel J (|^): 



27r(l-e-7*) 



exp 



-mf l {t)q_ l - mf 2 (t){q- f - q- t ) - m/ 3 (t) - <?_*) 

Here 7 = rj/m, q±i/f are initial/final values of q±, 



eT* - 1 

(15) 



f^t) = 1 fds Ids' I — coth ^ e -Ms-s') = 1 rt + 1 In 1 * (16) 
1 W 2 J J J 4tt 2T ' ' 2tt(T/uj c ) V ' 







and the functions f2(t) and fs(t) tend to the following values in the interesting limit of long times: 

mq 2 \ f dwwcoth^ 7 ^ oj c i T 



f 2 = (^-)= 7 ^^-^^l n ^: + _arctan-. (17) 

1 ' 2n u z + 7 Z 27r 7 7r 7 







and /3 = T + 2/2. The obtained exact solution allows to make several important observations. One of them is 
completely obvious: the particle q looses its coherence due to interaction with the Caldeira-Leggett bath of oscillators. 



Indeed in the long time limit we have f\(t) ^> f-x and the kernel (15) effectively reduces to 



J = l\I^Wf~ mhqlfS{q - l) " ^ mh9 - fS (9-i), (18) 

where L is the system size. In other words, for any initial conditions the density matrix tends to the same equilibrium 
form p(qi,q2) = (1/ ' L)e^ m ^ 2 ^ qi ~ q2 ' > which is not sensitive to the initial phase. According to eqs. (|ll|[l6|) the decay 
of off-diagonal elements of the initial density matrix is exponential at any nonzero T with the characteristic time 
To = l/riTq 2 ^ i . At T — the off-diagonal elements decay as a power law 

p i ( 9i )oc(^ c )-'N-* ) (19) 

but also in this case the information about the initial phase is practically lost in the long time limit. 

Another observation is that at sufficiently long times the average value of the kinetic energy of the particle mq 2 /2 
( |T7| ) is not zero even at T — irrespectively to its initial energy. At high temperatures T>7 ln(w c /7) the energy is 
given by its classical value T/2, but at lower T its value is determined by the interaction parameter 7 and the high 
frequency cutoff parameter lo c . It is sometimes believed that if initially all the bath oscillators are in the ground states 
and the particle energy is zero, no energy exchange between the particle and the oscillators will be possible because 
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the particle has no energy to excite the oscillators and the latter in turn cannot transfer their zero-point energy to 
the particle. This statement is obviously incorrect in the presence of interaction: the low energy particle will increase 
its average energy while the interaction energy will be lowered to preserve the energy conservation for the whole 
system. The presence of u c in eq. ([l7]) implies that there exists the energy exchange between the particle and all 
oscillators including the high frequency ones with uj ~ u c . One should not think, however, that such oscillators need 
to be excited in order to make this exchange possible. The energy transfer mechanism is different: in the presence of 
interaction the oscillator energy levels (including the ground state one) acquire a finite width and the oscillator can 
exchange energy in arbitrarily small portions. 

It is important to emphasize that some of the above effects cannot be correctly described within a naive perturbation 
theory in the interaction based e.g. on the Fermi golden rule. Just for an illustration let us choose the plane wave 
tp ~ exp(ipiq) as the initial state of the particle q and evaluate the transition probability W PlP2 to the state with the 
momentum p 2 . Without interaction one has W PlP2 = S PlP2 . For small 77 > at T — in the long time limit one gets 
from eqs. (J|), (|l5|) 

W P1P2 ~ — L^expf-- + (20) 

VxTx2 V v \xi xi) ) 

where xi — 41n(cj c £) and \2 — (2/ 71 ") -^(^c/t)- It is obvious that the result ((2^) cannot be recovered in any finite 
order of the perturbation theory in 77. 

The above model can also serve as an illustration of the role of "good" quasiparticles in the effect of quantum 
decoherence. It is clear that in this model one can carry out exact diagonalization of the Hamiltonian and introduce 
a new set of independent (and therefore coherent) particles/oscillators. It is also clear that this transformation will 
by no means influence our result for the density matrix p(q,q'): this result is exact. Thus also the calculation with 
"good" quasiparticles will yield the incoherent dynamics of q, however with much more efforts and with loss of physical 
transparency. The basic reason for dephasing of q is, however, transparent also in this case: q will be expressed as a 
sum of infinite number of independent particle/oscillator cordinates and therefore will be able to return to its initial 
state only after infinite time. 



C. Particle on a ring 

The second example is a quantum particle on a ring. Again we would like to calculate the probability W , but 
now we have to choose both qi and qf on a ring. We choose qi = qt = (see Fig.l), i.e. W is the probability for a 
particle to return to the same point q — 0. Again without interaction the two amplitudes K and K* decouple and 
can be evaluated separately. We have K = ^2 m K m , where K m ~ exp(iS m ) is the contribution from a path which 
traverses m times along the ring in a clockwise (m > 0) or a counterclockwise (m < 0) way and returns to the starting 
point q — 0. It is obvious that S m = S- m = 2ir 2 R 2 m 2 /t, where R is the ring radius. If we neglect terms KmK^, 
with ml 7^ ±m (one can argue that e.g. for S\ ^> 1 those are fast oscillating with m — ml terms which effectively 
cancel out in the course of summation over m and m'), then the probability W can be written as a sum of two terms 
W = W\ + W 2l where 

Wi = K ™ K * m > W 2 =J2 K m K*_ m . (21) 

m m 

The term W\ is determined by a pair of equivalent paths qi(t') — q2(t') (Fig. la). This contribution does not vanish 
in the classical limit. The term W2 comes from a pair of time reversed paths qi{t') = q 2 {t — t') (Fig. lb). This term 
describes the effect of quantum interference and therefore is very sensitive to the presence of the phase coherence in 
our system. Obviously, W 2 vanishes in the classical limit. As before, in both cases the phase factors exp(±iS' rra ) enter 
with opposite signs and cancel in each of the terms in (^l|). Without interaction we have W 2 — W\ — |-Ko| 2 - 

Now let us analyze the effect of interaction. The expression for the return probability W\ turns out to be insensitive 
to interaction due to exactly the same reason as in our first example: this probability is determined by the pairs of 
equivalent paths q\ = q 2 with F = 1. In contrast to Wi, the return probability W 2 determined by the pairs of time 
reversed paths is affected by interaction. We will restrict ourselves to the most interesting physical situation when 
the energy of a particle q remains conserved during its motion. Then the simplest pair of the time reversed paths 
is: q~(t') = 2i?sin(27rf'/£) and q+(t') is an even function of t'. Substituting these paths into ( |Io| , |ll"| ) after a simple 
integration we find 
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FIG. 1. A particle on a ring: (a) classical return paths, (b) time reversed return paths, (c) time reversed return paths of 
different sizes. 



2R 2 C 2 , / uj \ ( t Y sm\ut/2) 



The result (22) demonstrates that the return probability Wb for a particle </ after a time t for the time reversed paths 



acquires the factor 

F = exp(-Sj(t)), (23) 

due to interaction with a harmonic oscillator with a frequency w. The probability W\ remains unaffected. 

The whole consideration can be trivially generalized to the case of more complicated time reversed paths with an 
arbitrary winding numbers m. In this case we again find an additional factor ( |23| ) in the interference term W2 and 
Si is determined by an expression similar to $22\ ) which also depends on m. The term Sr is again zero for such paths 
for all m. Since no new effects emerge at to > 1, in what follows we will analyze only the simplest case m = 1 (p2]). 

The first conclusion one can draw from the result (^2|) is that no qualitative changes in the system behavior emerges 
if one varies the temperature T. The value Si(t) is smaller at low T <C uo as compared to the high temperature limit 
but the effect persists even at T — 0. Thus the relation between T and u> turns out to be important only in a 
quantitative sense, no qualitative dependence on this relation should be expected. 

A much more important parameter is Lot. We see that the value Si(t) and hence W2 oscillate in time with a period 
2ir/ui. After each such period the interference term W2 restores its "nonineracting" value while at all intermediate 
times the value W2 is smaller than in the noninteracting case. In this situation we still cannot speak about decoherence: 
the system keeps information about its initial phase and periodically returns to its initial state. On top of that if the 
radius of the ring R is constant in time (see below) we have Si cx l/t 2 , i.e. the oscillations of W2 practically disappear 
in the long time limit. All these results are not surprizing: one should not expect to find the decoherence effect in 
the system of two quantum mechanical particles. 

Quantum decoherence appears after the next step: we again modify our model coming from the interaction with 
one oscillator to the infinite set of oscillators with the ohmic spectrum ( |l3|) . In this case after the integration of (|22] ) 
over u> one gets 

riR 2 f sin 2 (7rz) ("kz\ , . 

Si = -!—\ zdz - L-i coth(— . (24) 



7T J (1-Z 2 ) 2 \TtJ 

This equation yields Si ~ r]R 2 at T = and Si ~ rjTR 2 t at Tt > 1. These results are in a nice qualitative agreement 
with those obtained in the exactly solvable model studied above. There the exponential decay of the off-diagonal 
elements of the initial density matrix p(q~i) oc exp(— t/ro) with To ~ was found at any T > while at T = 

a power law decay ( |l9|) was observed. In both cases the similarity is obvious if one interchanges R and q-i. 

We observe from (123) that if Si is small the influence functional F ~ 1 and no decoherence occurs. However for 
Si 3> 1 the interference of the time-reversed paths is completely suppressed, and quantum dephasing takes place. 
By setting Si ~ 1 we can, therefore, define the typical size of the ring R ~ L v beyond which the effect of quantum 
decoherence becomes important. In our particular example the characteristic dephasing length L v decreases as l/t 
at Tt ^> 1 and it is constant 
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L v ~ 1/y/rj (25) 

at T = 0. This is the effect of quantum decoherence. 

Note that by fixing R and increasing t in our problem we effectively decrease the particle velocity vq — 2nR/t 
which is eventually sent to zero as t approaches infinity. Therefore it is quite natural that at T = the interference 
term W2, although suppressed by a factor F — exp(— Si), does not decay in time. Qualitatively the same property is 
observed in our result for the decoherence time in a disordered metal (|l|): if we treat asa formal parameter which 
may be put equal to zero we will immediately arrive at a zero decoherence rate 1/t v in this limit. The same is true 
for the high temperature result [|l2|-|l5| . 

The above situation is, however, not very relevant for a metal where conducting electrons move with the velocity 
vq ~ vp which absolute value does not change in time. Thus to account for that we should rather keep the velocity 
of a particle i>o fixed. This implies that as we increase t the radius R for a classical return path increases linearly 
with time. In other words, we can slightly modify our model allowing our particle to choose the ring with a proper 
R for each time t (see Fig. lc). Then we immediately observe that Si ( p4| ) grows in time as Si cx t 2 at T = and 
Si oc t 3 at Tt ^> 1, i.e. in this case the decay of the quantum interference term W2 caused by interaction with the 
Caldeira-Leggett bath of oscillators is faster than exponential. Finally, let us note that in a disordered metal, although 
the length of the electron trajectory increases linearly with time, the dynamics is diffusive and the effective loop size 
grows as R 2 ~ Dt. Substituting this expression into (|5]) one obtains 

f exp(-t/r v ), t v = a/Dr), Tt < 1 

\ cx P (-(;/t v ) 2 ), T v =b(DrfT)- 1 ' 2 , Tt>l [ZG> 

where a and b are unimportant numerical coefficients of order one. 

Although the above simple model cannot be directly applied to disordered metals it demonstrates several important 
properties which will be also observed in a rigorous calculation. One such property is that the suppression of quantum 
interference between time reversed paths increases if the size of the loop grows in time. Another property is that 
oscillators with lj 3> T may give the dominating contribution to dephasing. In the above model the frequencies 
u> ~ 1/t give the maximum contribution, but if the spectrum of the problem is different (as e.g. in a disordered metal, 
see below) high frequency modes ~ lj c may also become important. The physical reasons for this conclusion were 
already clarified above: in order to have energy exchange in an interacting system it is not necessary to excite the high 
frequency oscillators, broadening of their ground state levels is sufficient. It is well known that the high frequency 
cutoff lo c enters the expression for the interaction induced decoherence rate of a quantum particle in the periodic [^]|| 
and the double well potentials jl0| . The same is observed here for an exactly solvable model of a free damped particle 
in the limit T = 0. 

Our model also demonstrates at which step of our calculation the effect of dephasing appears. Interaction of the 
particle q with one harmonic oscillator leads to the oscillations of the interference term Wz with the oscillator frequency 
u>. These oscillations are natural since the initial state is not the eigenstate of the system. Yet no dephasing appears. 
If coupling to many oscillators with different frequencies is introduced the probability W2 will be always suppressed 
and the system will never return to its initial state. Obviously this mechanism of dephasing has nothing to do with 
the temperature of the environment and it persists even at T = 0. In this limit the exponential decay of W% in time 



is due to increase of the loop size with t. The phase breaking length L v (25) depends only on the interaction strength 



and appears to be (roughly) insensitive to the particular (e.g. ballistic or diffusive) type of a particle motion. 



D. Scattering problem and perturbation theory 

One might wonder what is the relation between the above analysis and the standard perturbative treatment of a 
scattering problem for a quantum particle interacting with a harmonic oscillator. Beside its general importance the 
scattering approach is also of a practical relevance because it is frequently applied to conductance calculations in 
mesoscopic systems. 

The usual definition of a scattering problem operates with in- and out-scattering states measured after exactly 
infinite time. In other words, the limit t — * 00 is taken already at the very first step, and the whole calculation is 
carried out only in this limit. This is sufficient in many physical situations, but not for our problem due to the reasons 
to be clarified below. Here we will keep t finite (although possibly large) throughout the calculation and let it go to 
infinity in the end. With this in mind a direct connection to the scattering problem can be easily established. 

Consider a quantum particle q scattered on a harmonic oscillator with a frequency u>. Before scattering the oscillator 
is assumed to be in its ground state, and the particle is in the state with wave function ipi and the energy Ei. Assuming 
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the interaction to be of the same form as in (g) and proceeding perturbatively in the interaction strength C one can 
easily derive the probability W(t) for a particle q to leave the state ipi after the time t. It is just the sum of the 
transition probabilities Wu into all possible final states ijjf which are orthogonal to ipi in the absence of interaction. 
One easily finds (see e.g. Qj) 

W--£w,,.^\Wm? ^^m , (27) 

where AEif = u> + Ef — Ei. This result implies dW/dt oc sin(AEift) / AEft , i.e. at large t the transition rate 
experiences fast oscillations and approaches S(AEif) at t — > oo. If Et > Ei for all / (i.e. ipi describes the ground 
state of the noninteracting problem) the transitions are highly improbable at sufficiently long times and one may 
conclude that with the dominating probability the particle remains in its initial state and no quantum dephasing 
takes place. We would like to emphasize, however, that the time average of the escape probability W (|^) is not equal 
to zero already in this case. This is a direct consequence of the energy-time uncertainty principle. 

Let us now consider scattering on many oscillators. As before we will assume that the frequency spectrum for these 
oscillators is ohmic (|l3|). In order to define the scattering problem we also assume that interaction exists only in a 
certain space region —go/2 < q < ?o/2, i.e. we put C = Co(6(g + go/2) — Q(q — qo/2)). Without interaction the 
eigenstates of a problem are the plane waves ipi/f{q,t) = [L)^ 1 / 2 exp(ipi/fq — iE^/ft) where E^/f = p 2 ^/2m. The 
transition matrix elements can be easily evaluated. They give an important contribution for \pj — Pi\qo < 1 hi which 
case we get (i\Cq\f)\ 2 ~ Cq'ZoO 5 / ~Pi) 2 /L 2 . Subsituting this expression into ( |27| ) and carrying out the summation 
over the final states Yf ~~ > L J dpf an d over the oscillator frequencies (making use of ([l4|) with C — > Co), we find 

w® ~ 4. r ^ r 4 PM - P,f^p. 

L Jo J-i/qo \AE if y 

Further assuming that the initial particle energy is small Ei — > after simple integrations we obtain 

W(t) - ^ ]n((j c min(t, mq 2 )). (29) 

We observe that the escape probability W(t) is not zero and, moreover, it is not necessarily small. The reason for that 
is transparent. Although the contribution of each oscillator to W(t) is small, it is not zero at any finite t due to the 
energy-time uncertainty principle. The sum of these small contributions from many oscillators is finite and yields the 



result (29). We would like to emphasize that by no means this result is in contradiction with the energy conservation 
law, rather it demonstrates that one should be careful applying the energy arguments to describe the time evolution 
of an interacting quantum system, especially if it consists of an infinite number of degrees of freedom. Even at times 
much larger than the characteristic time scale (in our problem the relevant time scale is set by the dephasing time 
t v ) the energy uncertainty may be sufficient for W(t) to significantly differ from zero. 

In order to observe the relation of a perturbative expression (p9|) to our previous results we set q$ to be of order 



of the system size qo ~ L. After that the similarity between (|29|) and e.g. eqs. (19 2J) becomes completely obvious. 
Requiring that W(t, L) ~ 1 (the escape is complete) and using (|2^) with go ~ L one immediately arrives at the 
estimate ( ^5|) for the decoherence length L v ~ 1/ ^/rj derived previously with the quasiclassical analysis of the exact 
influence functional. This result leaves no room for doubts concerning the validity of the quasiclassical description 
of quantum dephasing at T = 0. In fact, the comparison of the results (|l9 29) with (|24|) demonstrates that the 



quasiclassical approximation rather underestimates the dephasing effect of interaction: L v may only become shorter 
if fluctuations around the classical trajectory are taken into account. 

Another obvious conclusion is that at any T (including T = 0) the effect of quantum decoherence in not only 
due to low frequency oscillators which, moreover, can be even completely unimportant. The relative contribution 
of oscillators with small uj depends on the particular form of the spectrum (~ w 7 ), being more important in the 
subohmic case 7 < 1 and practically irrelevant in the superohmic case 7 > 1 when the high frequency oscillators 
yield the main effect. The latter situation will be encountered in the next section where it will be shown that for a 
d-dimensional disordered metal one has 7 = 1 + d/2. Again we emphasize that the above results do not imply that the 
processes with high energy transfers are important for dephasing at any T. It is erroneous to interpret the parameter 
uj as describing the energy transfer: in our calculation the integral over u> always represents the summation over the 
bath oscillators (cf. eqs. (p3 - 14 )). No high frequency oscillators need to be excited, (small) energy uncertainty for 



infinitely many oscillators is sufficient to provide (large) dephasing for a particle q. 
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III. QUANTUM DECOHERENCE IN A DISORDERED METAL 



In order to provide a quantitative description of the effect of quantum decoherence in disordered metals it is 
necessary to go beyond the simple model considered in the previous section and account for (a) Fermi statistics and 
the Pauli principle, (6) the specifics of Coulomb interaction in a d-dimensional system and (c) the effect of disorder. 
The corresponding analysis is presented below. 



A. Density matrix and effective action 



Our starting point is the standard Hamiltonian for electrons in a disordered metal H c \ = Hq + H mt , where 

V 2 



H = J drV>+(r) 



2m 



H + U(t) 



Vv(r), 



(30) 



H int = \J dr J dr'V^(r)^(r')e 2 «(r-r')^(r')^(r). (31) 

Here p is the chemical potential, U(r) accounts for a random potential due to nonmagnetic impurities, and v(r) = l/|r| 
represents the Coulomb interaction between electrons. 
Let us define the electron Green-Keldysh function jl6| 

* = (Gu -G 12 \ _ [WTWiGv e iS KW 

\G 2l -G 22 ) /DWae^^l ' 1 ' 

where 

t 

iS^V^TrlnG^+i /dt' /dr ^l (Wz)2 . (33) 

o 

Here we performed a standard Hubbard-Stratonovich transformation introducing the path integral over a scalar 
potential field V in order to decouple the ^-interaction in ([H]) . In ( |32| , |33"| ) we explicitely defined the fields V\ (t) and 
V 2 (t) equal to V(t) respectively on the forward and backward parts on the Keldysh contour Jl(|. 
The matrix function Gy obeys the equation 

1 -KT - ^o(ri) + e^^i.ri)^ G v = 5{h - - r 2 ); (34) 



where 



H °- H0l -[0 ^-, + U(r))' V -{ V 2 (t,r))- ^ 
The solution of is fixed by the Dyson equation 

/ 

G v (hM) = G {h,t 2 ) - [ dt'G {h,t')eV(t')G v {t\t 2 ). (36) 



o 

The matrix Go is the electron Green-Keldysh function without the field V. 

It is well known that the 1,2-component of the Green-Keldysh matrix G is directly related to the exact electron 
density matrix 

p(t;r,r') = -iG 12 (t,t;r,r') = { Pv {t;r,r') Vl y„ (37) 

where we also defined the "density matrix" pv(t) related to the 1,2-component of the matrix Gv and performed the 
average over the fields V\ and V 2 as defined in ([52j). The density matrix p contains all necessary information about 
the system dynamics in the presence of interaction. 
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Making use of eqs. (|34|-|36|) after some formal manipulations (see M for details) one arrives at the equation describing 
the time evolution of the density matrix: 

i ~df = ~ eV+,pv ^ ~ C 1 ~ Pv)—^Pv ~ pv^j-O- - Pv), (38) 

where we defined V + = (Vi + V^)/2 and V~ = Vx — V2. It is important to emphasize that the derivation of this 
equation was performed without any approximation, i.e. the result (p8|) is exact. The equation (53) fully accounts for 
the Pauli principle which is important for the fluctuations of the field V~ . This field is irrelevant for dephasing. It is 
quite obvious from (Q) that the field V + (t, r) plays the same role as an external field. All electrons move collectively 
in this field, its presence is equivalent to local fluctuations of the Fermi energy fi — > fj, + eV + (t, r). The Pauli principle 
does not play any role here. There is no way how the density matrix p can "distinguish" the intrinsic fluctuating field 
V + from the stochastic external field, be it classical or quantum. Since the external field is known to lead to dephasing 
of the wave function, the field V + should produce the same effect. As the equation (p8|) is exact this conclusion is 
general and does not depend on approximations. 

In order to proceed further let us assume Coulomb interaction to be sufficiently weak and expand the action (|3J 
in powers of V up to terms proportional to V 2 . After a straightforward calculation (see e.g. W) one finds 



1 f d^fk kHme(u,,k) fur, 



2 J i2„)* V ~(-">- k) t ' ' <"M-jy-(a,,fc); (39) 
where e(u>,k) is the dielectric susceptibility of a disordered metal 

e(u,k) = l+ . 4 ™ 2 - (40) 
— ltd + Dk £ 

Here a = 2c 2 NqD is the classical Drude conductivity, Nq is the metallic density of states and D = uj?Z/3 is the diffusion 
coefficient. For the sake of simplicity in eq. (|4(]) we disregarded the phonon contribution which will not be important 
for us here. The expression ( ^0| ) is valid for small wave vectors k < l/l and small frequencies uj < l/r e = vp/l. 

Note, that if one considers only nearly uniform in space (k w 0) fluctuations of the field V one immediately observes 
that eqs. (^9|,^0|) exactly coincide with the real time version of the Caldeira-Leggett action |5|,^|[| in this limit (cf. 
(p| |n| , |l4| )). Taking into account only uniform fluctuations of the electric field one can also derive the Caldeira-Leggett 
action expressed in terms of the electron coordinate only. In this case the effective viscosity 77 in the Caldeira-Leggett 
influence functional is 77 ~ e 2 R s , where R s is the sample resistance (in contrast to the effective viscosity for the field 

V which is proportional to 1/R S ). 

For our present purposes it is not sufficient to restrict ourselves to uniform fluctuations of the collective coordinate 

V of the electron environment. The task at hand is to evaluate the kernel of the operator 

J = Y J U\V){V\U + 

v 

where the sum runs over the states of the electromagnetic environment with all possible k and ui. Averaging over 
V + , V~ amounts to calculating Gaussian path integrals with the action (|3^) and can be easily performed. We obtain 

ri(t)=r 1/ r 2 (t)=r 2/ 

J(t,t';r 1 f,r 2 f 1 r 1 i,r2i) = J Dn J Dr 2 / ^DpiDp 2 x 

ri(t')=rij r 2 (t')=r 2 i 

x exp {iSb[ri,pi] - i5 [r 2 ,p 2 ] - iS R [ri, pi, r 2 , p 2 ] - 5/[ri,r 2 ]}; (41) 

where 

t 

(42) 



S [r,p] = J dt"Uv-^--U{v) 



is the electron action, 
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t t 

,2 



Sfl[ri,pi,r2,p2] = ^ J dii / d* 2 {JZ(ti - t 2 ,ri(*i) - n(* 2 ))[l - 2n(pi(t 2 ),ri(i a ))] - 
t' t' 

-iZ(ix - i a , r a (ti) - r 2 (t 2 )) [l - 2n(p 2 (i 2 ), r 2 (i 2 ))] 
+i2(ti - ia, ri(ix) - r 2 (i 2 )) [l - 2n(p 2 (i 2 ), r 2 (t 2 ))] - 

- J R(ii-i 2l r 2 (ii)-r 1 (i 2 ))[l-2n(pi(i 2 ),r 1 (i 2 ))]}, (43) 
where n(p, r) is the occupation number and 

2 * * 

5j[r l!r2 ] = y y dii y diaj/^x - <a,n(*i) -ri(*a)) + /(ti- i 2 ,r 2 (ii) - r 2 (i 2 )) - 

-/(ii - t a ,r!(ti) - r 2 (i 2 )) - /fa - t a ,r 2 (ti) - r x (i 2 ))|. (44) 

In equilibrim n is just the Fermi function. In this case at the scales |r| > I the functions R and I are defined by the 
equations 

'0, r) = / ^?lm ( co«h ( ^ e -«+«" (46) 



(2tt) 4 \k 2 e{uj,k)J \2T) 

The expression in the exponent of eq. ([Il]) defines the real time effective action of the electron propagating in a 
disordered metal and interacting with other electrons. The first two terms represent the electron action Sq ( f42| ) on 
two branches of the Keldysh contour while the last two terms Sr and Si determine the influence functional (cf. eqs. 
©ErO'l^) °f th° effective electron environment. As can be seen from eqs. (j43] - ^6|) this influence functional is not 
identical to one derived in the Caldeira-Leggett model. However on a qualitative level the similarity is obvious: in 
both models the influence functionals describe the effect of a certain effective dissipative environment. 



B. Decoherence time 



Let us first neglect the terms Sr and Sj describing the effect of Coulomb interaction. Then in the quasiclassical 
limit ppl ^> 1 the path integral ( f4l| ) is dominated by the saddle point trajectories for the action Sq: 

p = -dHo/dr, f = dH Q /dp (47) 

with obvious boundary conditions ri(i') = r^, r(t) = rif for the action 5o[ri,pi] and r 2 (t') = r 2 j, r 2 (i) = r 2 / for 
the action S'o[r 2 ,p 2 ]. 

Since in a random potential U(r) there is in general no correlation between different classical paths ri(t) and r 2 (i) 
these paths give no contribution to the integral (pl|): the difference of two actions So in the exponent may have an 
arbitrary value and the result averages out after summation. Thus only the paths with So[ri, Pi] — 5 , o[r 2 , p 2 ] provide 
a nonzero contribution to the path integral (]4l|). Two different classes of paths can be distinguished (see e.g. 

i) The two classical paths are almost the same: ri(i") ~ r 2 (i"), Pi(i") — p 2 (i") (cf. Fig. la). For such pairs we 
obviously have ru ~ r 2i and ri/ ~ r 2 j . Physically this corresponds to the picture of electrons propagating as nearly 
classical particles. In the diffusive limit these paths give rize to diffusons (see e.g. | jl2l ) and yield the standard Drude 
conductance. 

ii) The pairs of time reversed paths. In this case — r 2 /, r%f ~ r 2 ^ (cf. Fig. lb,c). In the path integral (||) the 
trajectories ri and r 2 are related as r 2 (i") ~ rx(i + 1' — t") and p 2 (i") — — pi(i + 1' — t"). In a disordered metal these 
paths correspond to Cooperons and give rize to the weak localization correction to conductivity Sad- This correction 
is expressed in terms of the time integrated probability W(t) for all diffusive paths to return to the same point after 
the time t (see e.g. [|l2| , ^3| ). In the absence of any kind of interaction which breaks the time reversal symmetry this 
value coincides with the classical return probability and is given by the formula W (t) = (4nDty d / 2 a~^~ d \ where d 
is the system dimension and a is the transversal sample size. 
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The weak localization correction Sad diverges for d < 2. This divergence can be cured by introducing the upper 
limit cutoff at a certain time t v . This time is usually reffered to as decoherence time. One finds [ jijLjill : 

s<r d = { y u;> (48) 

The physical reason for the existence of a finite t v is the electron-electron interaction which breaks the time reversal 
symmetry in our problem. To evaluate t v we first note that the functions R and I ( [i"6|) change slowly at distances of 
the order of the Fermi wavelength 1/pf- Therefore we may put ri(i") = r(t"), i"2(t") = Yt{t") = r(t + t' — t"). Here 
r(t") is a classical trajectory with the initial point r(t') = and the final point |r(£)| < I, i.e. we consider trajectories 
which return to the vicinity of the initial point. The contribution from the time reversed paths Wi to the return 
probability has the form 

W 2 (t - t') ~ W 20 (t - t')F, F = ^ e -iSn[t,t';r,p;r i .,p i .]-S I [t,f;r,r i .]^ ^ ^ 

where W20 is the probability without interaction and the average is taken over all diffusive paths returning to the 
initial point. The value F in (^) decays exponentially in time, therefore we may put the average inside the exponent. 
It is easy to observe that the term Sr gives no contribution to this average. Working out the average of Si we obtain 

+00 

F = e -<S/[t,0;r,r«]) r = exp I _ fe 2 f d ^/(^ r ( t j _ r ( ))^ I ( 50 ) 



To find the average over the diffusive paths, we introduce the Fourier transform of the function I(t, r) and replace 

( e -ik(r(t)-r(t'))) r by e -Dfe 2 |t-t\ Then we get 

te 2 f dud d k ( -4vr \ , / a \ Dk 2 , s 

lni? (*) = --*-7 / 7^TT Im tt-t-it coth — - 2 ■ n2> . (51) 



(2 7 r) d + 1 \k 2 e{uj 7 k)J \2T ) uo 2 + D 2 k A 

Let us first consider a quasi-one-dimensional system with a < I. Making use of eqs. (|5l],^) and integrating over k 
we find 

lnF(*)=-t^ y^coth(./2T)_ 

(71 J 27T ^/W 

1/t 

The upper cutoff in ( p2| ) is chosen at the scale ~ l/r e because at higher u> the diffusion approximation becomes 
incorrect. From (52) we obtain 

Jcxp(-i/%), T v ^7r{3/2) 1 / 2 a 1 /e 2 v F , t « 1/T 2 r e 

* ~ \ exp(-(i/r v ) 3 / 2 ), r v ~ 0.5(^ 1 /Z? 1 /2 T) 2/3 ] 4 » 1/T 2 Te l^J 

In the long time limit (which cannot be reached at T — ► 0) the decay of F is faster than exponential (cf. p5fl). At 
sufficiently low T < (rgT^,) -1 / 2 this difference becomes important only at t r v where F is already exponentially 
suppressed. For smaller i the decay is exponential and the decoherence rate increases linearly with T. Defining the 
effective decoherence length L v ~ ^/Dt^, in the temperature interval 1/2t v (0) < T < (leT^O)) -1 ' 2 one finds 

2 V2^4^T, (54) 



L 2 (T) L 2 (0) TrDcn 

where i- v (0) is the dephasing length at T = 0. At temperatures lower than l/2r v (0) the decoherence time and 
length are almost temperature independent. With the aid of eqs. ( |4^ , |5^ ) it is also easy to find the weak localization 
correction Sa\ to the Drude conductance. In the limit T — we obtain 

^ = -^- * (55) 
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i.e. Sai ~ — (Ti/^/TVch, where iV c h ~ PpS is the effective number of conducting channels in a Id mesoscopic system and 
<7d = <ra 3 ~ d is the Drude conductance of a <i-dimensional sample. Thus for iV c h 3> 1 the weak localization correction 
is parametrically smaller than the Drude conductance even at T = 0. 
For 2d and 3d systems the same analysis yields 



1_ _ e 2 

t v 4:TTa 2 T e 



[l + 27Y e ln(TTv)], 2d, 



- N3/2 



\ + &{TT e yi% 3d, (56) 



The influence functional F decays exponentially F = exp(— t/r^) except for the case of 2d systems at high T where one 
has ln_F cx —tint. We observe that also in this case the difference from a purely exponential decay is not significant 
and can be ignored. 



C. Further remarks 



The above formalism provides a rigorous description of the effect of quantum decoherence in disordered metals and 
allows to derive the decoherence rate 1/t v at low T for various dimensions. The physical reasons for dephasing of 
electron wave functions remain the same as in the case of a particle interacting with the Caldeira-Leggett bath of 
oscillators. The corresponding discussion is presented in the previous section. Here we only add several comments. 

It might be interesting to investigate the average kinetic energy of an electron in a disordered metal in the presence 
of interaction. At low T it turns out to be temperature independent due to the same reason as for a quantum particle 
in the bath of oscillators, i.e. due to interaction. A naive calculation along the same lines as in Section 2 yields 
(mr 2 /2) — (j, ~ 1/t v at T = 0. A more accurate analysis is beyond the frames of the present paper. This result is just 
a manifestation of the well known fact: the distribution function n(p) of interacting electrons (not quasiparticles) is 
nonzero for any p even at T = (ItJ . 

An obvious consequence of our results is that the decoherence rate l/r v exceeds temperature at sufficiently low T. 
Does this fact imply the breakdown of the Fermi liquid theory (FLT) hypothesis at such Tl From a formal point 
of view the violation of the inequality Tt v > 1 is yet insufficient for such a conclusion (see also jl8|), although it 
obviously does not support the FLT hypothesis either. It was demonstrated in Section 2 that a finite decoherence 
rate for a certain variable of interest has no direct relation to the lifetime of quasiparticles (the latter is infinite in 
the case of a free particle in the Caldeira-Leggett bath) . In our calculation t v determines the dephasing time for real 
electrons and not for the Landau quasiparticles. Our result implies that interacting electrons are "bad" particles since 
their wave functions dephase even at T = 0. The possibility to construct "better behaving" quasiparticles remains 
questionable for disordered metals. Furthermore, even if such quasiparticles exist their properties are completely 
unknown. Therefore at this stage they can hardly be used for calculation of any physical quantity. An important 
advantage of our method is that it allows for a direct calculation of measurable quantities in terms of interacting 
electrons without appealing to the Fermi liquid hypothesis. 



IV. DISCUSSION OF EXPERIMENTS 



Our results for the decoherence time at low T turn out to be in a remarkably good agreement with available 
experimental data obtained in various physical systems. In Ref. p| we have already carried out a detailed comparison 
between our theory and the experimental results jlj obtained for Id Au wires. The agreement within a numerical 
factor of order one was observed for all samples studied in [fij . Since both t v and L v are defined with such an accuracy 
no better agreement can be expected in principle. 

Here we will present a comparison of our theory with two other experiments carried out with 2d electron gas in 
semiconductor structures. These systems are somewhat different from the metallic ones mainly because of much 
higher effective resistance. The parameters of the systems were chosen in a way to realize a quasi Id conducting 
system with disorder. One such experiment was carried out by Pooke et al. fl9|| in narrow Si pinched accumulation 
layer MOSFETs. These authors studied samples with resistances ranging from 120 to 360 kSl and observed a finite 
decoherence rate 1/t v at all temperatures. On a log-plot (not shown) clear signs of saturation at low T are seen. 
The corresponding data for L v [fill obtained for 3 samples are presented in Fig. 2 and in the Table 1 together with 
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our theoretical predictions. Note, that the experimental value of L v at T = was obtained by a linear extrapolation 
l/L 2 9 {T) = 1/L*(0) + BT of the experimental data @. 

Table 1 



Vg, V p , V 


L-p(0), ^m 


L^ COT (0), M m 


B cx p, t im- 2 K- 1 


B thcm \ fim-' z K- L 


60, 


1.24 


0.78 


5.45 


3.8 


60, -3 


0.85 


0.58 


7.83 


5.4 


75, -6 


0.66 


0.53 


6.65 


6 




FIG. 2. The experimental data (symbols) fitted to our theory (solid lines). 

The agreement between theory and experiment is within a numerical factor of order one, i.e. again within the 
accuracy of the definition of L v . Adjusting this numerical factor we observe a perfect fit of the experimental data by 
our theoretical curves (see Fig. 2). As predicted (cf. eq. (54)) at sufficiently low T the value l/L^ increases linearly 
with temperature. 

Very recently new measurements of the dephasing time t v in quasi Id <5-doped GaAs structures with few conducting 
channels were reported pO]j . The typical effective resistance of the samples was 6-h 30 Mfi, i.e. it was even higher 
than in Ref. [l9| , the length and the width of the wires were 500 /j,m and 0.05 /im respectively. Experimental results 
for 3 different values of the gate voltage V g (20) are presented in the Table 2 together with our theoretical predictions. 
Again a perfect agreement between the maximum measured value L v and the one derived from our theory is found 
for all V g . 

Table 2 



Vg,V 


n, 10 12 cm- 2 


r, Mn 


L%*, Atm 


L^ oor (0), f im 


+0.7 


4 


5.94 


0.35 


0.23 





2.7 


~ 18 


0.09 


0.08 


-0.35 


2 


~ 36 


0.06 


0.04 



The temperature dependence of the data L v |2Cj is also in excellent agreement with our predictions at all T, see 
Fig. 3. 

In Ref. po| the measured maximum values for t v were compared with the theoretical formula suggested in Refs. 
ID, pl| and a discrepancy by a factor of 50 was reported. This fact allowed the authors to conclude (20j that their 
experimental results for argue against the idea of decoherence by zero-point fluctuations of the electrons. The 
comparison of the data pQ] with our theoretical results clearly demonstrates that this conclusion is simply an artefact 
of the inadequate choice of a theoretical formula adopted in EQ]. The measurements of the decoherence rate M| 
strongly support the idea of decoherence due to intrinsic quantum noise rather than argue against it. 

In Ref. [M also another interesting experimental observation was made: a crossover to a highly resistive state 
was found around T~ 1 K. This observation was interpreted in p0| as a Thouless crossover to the regime of strong 
localization and was also qualified as contradicting the very idea of quantum decoherence at T = 0. Several comments 
are in order. 

(i) For the systems studied in |^| one has ^/N^ ~ 2 -j- 3. Hence, the weak localization correction at low T ( |55| ) 
should be of the same order as the Drude conductance (in fact, this is exactly what was observed in p0|]) and the 
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1 10 

T, K 

FIG. 3. The experimental data Q (symbols) fitted to our theory (solid lines). 

Thouless crossover cannot be ruled out theoretically for such small N c ^. Also the result for t v may become more 
complicated for small JV c h [0 because of a somewhat more important role of capacitive effects. This may in principle 
lead to a relative increase of L v (it appears, however, that this effect is not very important in |pof ). Thus even the 
presence of a Thouless crossover in the samples [|o| would by no means contradict our theoretical results and hence 
the idea of quantum decoherence at T = 0. 

(ii) The interpretation of the observed effect as a Thouless crossover is not quite convincing. This interpretation 
is based on two reasons p0fl : (i) at the crossover L v was found to be "only" ~ 3 times smaller than Zi oc and (ii) the 
resistance of a wire segment of the length ~ Zi oc was found to be ~ 20 kfi. The first observation does not contradict 
to our theory which predicts l\ QC jL v ~ V-^ch ~ 3 at low T. As to the resistance of a segment ~ Zi oc it is of the order 
of the quantum resistance unit also at temperatures T ~ 3 -j- 10K, i.e. well above the crossover. Hence, no definite 
conclusion can be drawn. 

(iii) The expression for l\ oc used in |m| makes sense only provided the condition L v > L li oc is satisfied (L is the 
wire length). For L v <C L the wire can be viewed as N ~ L/L v independent samples connected in series (in |2(J N is 
typically of order 10 3 ~ 10 4 ) . Since the the length of each of such samples L v is ~ 3 times smaller than the localization 
length li oc it is somewhat naive to seriously believe that strong localization can be observed in such samples relying 
only on the "order-of-magnitude" character of the relation between L v and Z/ oc at the crossover. Moreover, assuming 
the dependence L v oc T- 1 / 3 one immediately observes that the "standard" crossover condition Zi oc ~ L v would hold 
at temperatures ~ 30 times smaller than the actual crossover temperature, i.e. at T ~ 20 -j- 30 mK. We see no way 
how the Thouless crossover can be expected above this temperature range for the systems studied in [f20"f . 

(iv) The inequality Tt v ^ 1 is violated in |2(]] at all relevant T: at the crossover one has Tt v ~ 0.3-^0.5 depending 
on the sample, and Tt v remains of order one even at T ~ 10 4- 20K, i.e. well in the weak localization regime. If 
one believes that the violation of the above condition signals the breakdown of FLT due to interaction the whole 
discussion of the Thouless crossover becomes pointless. An alternative would be to acknowledge that the value Tt v 
may not be the relevant parameter as far as FLT is concerned. But in any case real electrons dephase and therefore 
can hardly be described within the Thouless scenario of strong localization. 

(v) Since Tt v < 1 and Zi oc is several times larger than L v at the crossover interaction definitely plays a very 
important role "helping" to localize electrons instead of destroying localization. If so, why not to assume that the 
whole effect is solely due to interaction and not due to spacial disorder? For instance, it is well known that the 
mobility of a quantum particle in a periodic potential can decrease dramatically with T if this particle is coupled to a 
dissipative environment (see e.g. ||). At T = this particle can even get localized due to the effect of quantum noise 
of the environment |^,|| which completely destroys the phase coherence of the wave function. Within this scenario 
the crossover to a highly resistive state p0[ can be considered as supporting the idea of quantum decoherence due 
to intrinsic quantum noise. Thus, although a detailed interpretation of the crossover ]2C| ] is still an open problem, 
presently we see no way to use it as an argument against quantum decoherence at T — 0. 

It was argued in p0| that the effect of saturation of t v observed in many experiments at low T can be caused by 
the external microwave noise. Although filtering of external noise is indeed a serious experimental problem it is quite 
obvious that the above explanation faces several severe problems. Without going into details let us just indicate some 
of them. 
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Firstly, according to the arguments |20j the dephasing effect of the external noise may not be accompanied by 
heating for low resistive samples with R <C 24 kfi, whereas in the opposite case of highly resistive samples heating 
is unavoidable. It is not clear how to match this conclusion with the experimental results jl9| where clear signs of 
saturation of L v at T < IK were seen for samples with resistances up to R ~ 360 kSl ^ 24 kil. 

Secondly, the formula (Q) quantitatively (within a factor of order one) describes the low temperature value of the 
decoherence time t v measured in different experiments in at least 10 Id semiconductor and metallic samples. In order 
to interpret the results of all these measurements in terms of external noise one should assume that external noise 
always adjusts itself to a particular value of the sample conductance (in various experiments these values differ by 
many orders of magnitude) and the Fermi velocity. More than that, the corresponding electric field produced by the 
external noise inside the sample should be always of the same order as one due to the intrinsic quantum noise. It 
would be interesting to estimate the probability for such a coincidence in (at least) 10 different samples. 

Thirdly, the presence of the external noise can only be proven experimentally by making experiments with and 
without necessary filtering and observing different results in these two cases. In Ref. |2p| the external noise power 
was estimated, however no evidence for its existence in experiments was presented. In contrast, quantum noise with 
uj > T is well observable reality, see e.g. 

Thus the explanation of the existing experimental results in terms of external noise turns out to be problematic. 
Perhaps new experiments are needed to unambiguously rule this issue out. 

V. CONCLUSIONS 

In the present paper we have discussed the fundamental effect of interaction induced quantum decoherence in 
disordered metals. 

We have considered a simple model of a quantum particle interacting with the Caldeira-Leggett bath of oscillators. 
An exact solution in the case of a free damped particle demonstrates that the off-diagonal elements of the particle 
density matrix decay in the long time limit at all temperatures including T = 0. For a particle on a ring similar 
results are found. 

A very transparent physical picture of the effect of quantum decoherence due to interaction with a quantum bath 
of oscillators emerges from our analysis. The interference contribution to the return probability Wi for a particle 
interacting with one oscillator with a frequency ui oscillates in time and is smaller than one for all time moments 
except t — 2ixnju} when the system returns to its initial state. If interaction occurs with infinitely many oscillators 
with a continuous distribution of frequencies the particle will never return exactly to its initial state. At T = the 
return probability will be suppressed by a factor 

W 2 cx exp(-ni? 2 ), 

where rj is the viscosity of the environment and R is the size of the return path. This defines the typical dephasing 
length L v ~ l/^/rj. For a particle in a diffusive environment the typical size of such a path grows with time as 
R ~ y/Dt, and the interference probability will decay as W2 oc exp(—r]Dt). This is the effect of quantum dephasing. 
It has an essentially quantum mechanical nature, therefore its existence at T — is by no means surprizing. At 
nonzero T the effect increases due to increasing fluctuations in the environment. 

The main features of the effect derived by means of our simple model are reproduced within the rigorous analysis 
developed for a disordered metal, in this case one should account for the Pauli principle and the sample dimension. 
The electron interacts with the fluctuating electromagnetic field produced by other electrons. This field can be again 
represented as a collection of oscillators with a somewhat more complicated spectrum than in the Caldeira-Leggett 
model. Quantitatively the results will depend on that, but the physical nature of the effect remains the same. 

The effect of quantum dephasing studied here has no direct relation to the question about the existence of "coherent" 
quasiparticles in the problem. It is obvious that only the behavior of measurable quantities is of physical importance. 
If these quantities (e.g. the current or conductance in the case of disordered metals) are expressed in terms of 
"incoherent" variables, quantum dephasing yields directly measurable consequencies. A convincing illustration for 
that is provided by the existing experimental data. 

Our results indicate a necessity to reconsider the commonly adopted point of view on the role of interactions in 
disordered metallic systems. In particular, one arrives at the conclusion that no electron localization takes place and 
low dimensional disordered metals with generic parameters do not become insulators even at T = 0. 
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